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Abstract 



We report on numerical studies of the influence of Coulomb interactions 
on localization of electronic wavefunctions in a strong magnetic field. Interac- 
tions are treated in the Hartree-Fock approximation. Localization properties 
are studied both by evaluating participation ratios of Hartree-Fock eigenfunc- 
tions and by studying the boundary-condition dependence of Hartree-Fock 
eigenvalues. We find that localization properties are independent of interac- 
tions. Typical energy level spacings near the Fermi level and the sensitivity 
of those energy levels to boundary condition show similar large enhancements 
so that the Thouless numbers of the Hartree-Fock eigenvalues are similar to 

those of non-interacting electrons. 

PACS numbers: 71.30. -hh, 71.55.Jv, 73.40 Hm. 
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In studies of the integer quantum Hall effect it is usually assumed without justification 
that Coulomb interactions between electrons can safely be ignored. At strong magnetic 
fields, eigenfunctions of the single-particle Hamiltonian are localized by a disorder potential 
at almost all energies and the localization length diverges as |-E — Ecn]'"^^^ at a discrete set 
of critical energies E^n near the center of each disorder-broadened Landau level []I|. The 
critical energies are similar to the mobility edges of the Anderson metal-insulator transition, 
except that the localization length is finite both for E > Ecn and for E < Ecn- Theory 
suggests that at T = the Hall conductivity {(Txy) jumps by e'^/h each time the Fermi energy 
Ep crosses one of the critical energies, that the dissipative conductivity {crxx) is zero unless 
Ep = EcN, and that a^x = (^xy = = when Ep = E^n- All these suggestions are 
supported by experimental studies, except possibly 0] for the predicted universal value of 
the critcal dissipative conductivity axx- In this Letter, we report on numerical studies of the 
integer quantum Hall effect which include Coulomb interactions at the lowest consistent level 
by performing Hartree-Fock calculations in the presence of random disorder potentials. The 
possible importance of Coulomb interactions is suggested by studies in three-dimensions 
at zero magnetic field where the metal-insulator transition can be profoundly altered by 
interactions and by the presence of a Coulomb gap in the integer quantum Hall regime 
PJlO(| . However, general scaling considerations suggest that for the quantum Hall effect 



the critical conductivities are independent of Landau level index A^, details of the disorder 
potential, and of interactions. 

The concept of a localization length can be generalized to the case of interacting elec- 
trons by defining it in terms of the dependence of the imaginary part of the the one-particle 
Greens function on spatial coordinates. In the Hartree-Fock approximation, this change 
is accomplished by replacing the single-particle eigenfunctions in a random potential by 
the self-consistent Hartree-Fock eigenfunctions. We have studied the localization proper- 
ties of these eigenfunctions by examining their participation ratios and find no qualitative 
changes. This result is reasonable since the self-consistent Hartree-Fock Hamiltonian is 
always equivalent to the single-particle Hamiltonian for a non-interacting system with an 



effective random potential which however presumably has smaller long wavelength compo- 
nents because of screening effects. (The exchange potential is non-local but its effect is 
equivalent [|1T| to that of a local potential.) Universality of localization properties in the 
absence of interactions would then seem to imply no change in critical localization properties 
in the Hartree-Fock approximation. 

The influence of interactions on the critical conductivity is difficult to determine, even 
within the Hartree-Fock approximation. The numerical evaluation of axx necessitates the 



inclusion of interaction vertex corrections by summing ladder diagrams [|12| in the presence 
of disorder. We have avoided this considerable obstacle by appealing to the close relationship 
between a^x and the more easily calculated Thouless numbers which are defined in terms of 
the dependence of eigenenergies on the boundary conditions of the system. The definition 
we use for the Thouless number is: 

where 1/AE = L'^D{E) is the average level spacing calculated from the density of states per 
unit area D{E), {6E) is the geometric mean of the eigenvalue differences between periodic 
and antiperiodic boundary conditions and is the area of the finite system. For non- 
interacting electrons the Kubo formula for the conductivity can be written as 

axx = ne'LW{Ep){\Vx\') (2) 



where is an average |]T3[ over matrix elements of the velocity operator between states 

at energies near Ep. It follows from simplifying assumptions about correlations between 
velocity matrix elements and energy level spacings that 

..^^^ (3) 



both at [|T3| B = and at p!^p!5| strong fields. (In the metallic limit a rigorous relationship 



between axx and a suitably defined Thouless number can be established at i? = 
for non- interacting electrons.) Our numerical calculation is motivated by the expectation 



that within the Hartree-Fock approximation a quahtative relationship between QTiEp) and 
a^x will survive interactions. Hartree-Fock eigenvalues change with boundary conditions 
p!3| , p^ directly because of the velocity matrix elements between Hartree-Fock eigenfunctions 
and indirectly because of the change in the exchange potential produced by the changed 
eigenfunctions. The former effect is related to the 'bubble' diagram for the conductivity just 
as it is without interactions, while the latter effect is related to vertex corrections in the 
Hartree-Fock approximation for the conductivity. In the integer quantum Hall regime D{Ef) 
is strongly suppressed by interactions [^|lOl just as it is for disordered electrons at 5 = 0. 
In view of Eq.( |^), this suggests that o"c could be different in interacting electron systems. 
We have found that a remarkable cancellation occurs in which the change in gxiEp) due to 
the decrease in D{Ep) (increase in AE) is cancelled by increased sensitivity of the Hartree- 
Fock eigenvalues to boundary conditions. To the precision of our calculations, gx^Ep) is 
unchanged by interactions. 

In the model used for our numerical calculations N electrons in the lowest Landau level 
are confined to a square of area = 27r£^A^^ = N^^q/B. is the number of single- 
particle states in the lowest Landau level, $0 = hc/e is the magnetic flux quantum and the 
magnetic length, C. is defined by these relations. Our model disorder consists of a randomly 
located delta-function 'impurities' with a random strength uniformly distributed between 



—A and A. For this model [TJ] the energy scale which characterises the Landau level width 
is r = (X^Nj/PL'^Y^'^ where Nj is the number of impurities. The relative strength of 
Coulomb interactions and disorder is specified by the parameter 7 = {e^/eVj/T {e is the 
dielectric constant). The solution of the Hartree-Fock equations is greatly facilitated by the 



fact that matrix elements of the exchange potential [TT|] can be expressed in terms of the 
electron density rather than the density matrix; technical details of these 'self-consistent 
field' calculations have been explained elsewhere. p!7|,p| 

Localization properties were investigated by evaluating participation ratios for self- 
consistent Hartree-Fock eigenfunctions, (pa- 
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^ [/^r>.(r)|f 

Figure |I| displays InP^ as a function of z/ for 7 = and 0.4 and up = 1/2. (We reserve the 
symbol for the filling factor at the Fermi level while u denotes the fraction of the density 
of states below a particular energy. Most of the results reported here were calculated for the 
case Up = 1/2 for which the extended state occurs at the Fermi energy.) The participation 
numbers are shifted slightly upward by interactions near u = 1/2 and appear to have a cusp 
in their dependence on u (Note that no cusp appears in the dependence on energy since 
the linearly vanishing density of states implies that |z/ — l/2\ oc [E — E^)"^ in the presence 
of interactions.) The fact that participation numbers are not greatly changed when the 
localization length is smaller than our finite system sizes suggests that localization lengths 
are quantitatively unchanged by interactions. We infer from this that electron-electron 
interactions do not affect the critical properties of localization in any significant way and 
that the exponent remains 7/3, in agreement with recent experiments 0. 

To examine the behavior of the localization length in more detail we have investigated 
the system-size dependence of our results. The left panel of Fig. 2 shows the dependence of 
InPa on InA^^ in the absence of electron-electron interactions for z/ = 0.5, z/ = 0.066 and 
V = 0.0022. The right panel of Fig.2 shows similair results in the presence of electron- 
electron interactions. Results are calculated at the filling factors 0.5, 0.081 and 0.056 for 
Up = 1/2 and at 0.538,0.147 and 0.019 for up = 1/5. For localized states we expect that 
P oc {i/L^ oc S,"^ /N^. The results displayed in the left and right panels show that states 
in the Landau level tails are well localized whether the Fermi energy is at E^. or located in 
the Landau level tail. The results shown at z/ = 0.5 and 0.538 are in a regime where the 
localization length exceeds the system size. In this regime the participation ratio for non- 



interacting electrons is known ||T8[ to have a powerlaw dependence on system size, P oc A^^ ^ 



where the exponent A is related to the multifractal character |jT9[ of the extended states 



responsible for the anomolous diffusion which occurs when Ep = E^. The results in 
Figure show that the self-consistent Hartree-Fock eigenfunctions in the critical regime have 



the same fractal properties as the critical wavefunctions for non-interacting electrons. A is 
given by D[2\ = 2(1 — A) where D[2] is a measure of the structure of the multifractal. From 



Fig. 1^ we estimate that D[2] ~ 1.6 in agreement with earlier work [|I8],|T9|. 

In Figure ^ we show, T/AE, 5E/T and qt as a function of v for vp = 1/2, A^;^ = 72 
and Nd = 1412 for both interacting and non-interacting electrons. Our results for non- 
interacting electrons are in agreement with earlier [|14|,|15] work. As expected, interactions 
lead to an increase in AE near Ep corresponding to the density-of-states suppression. How- 
ever, 6E is also increased and the net result is that the Thouless number is unchanged to 
within the accuracy of our calculations. In agreement with Hanna et al. we find that the 
peak Thouless numbers appear to already approach their thermodynamic limit for A^^ ~ 50. 
Our numerical results suggest that the peak Thouless numbers are independent of interac- 
tions. The Thouless numbers away from Ep are larger for finite-size systems in the presence 
of the interactions, presumably because the screening changes the nature of the disorder po- 
tential. For bigger system sizes we expect that AE oc 1/L at Ep which implies that 6E 
is also (X 1/L. If Eq.( ^ were exact, the critical conductivity would be proportional to the 
peak Thouless number. Our Thouless numbers, both with and without interactions, would 
correspond to peak conductivities ~ 0.2(e^//i). This value is smaller than both the expected 
critical conductivity /2h and the self-consistent Born approximation peak conductivity, 
e^/ZiTT. To understand these quantitative relationships the connection between gp and a^x 



must be established under conditions more general that those considered in . 

In conclusion, we have shown that for the metal-insulator transition of the integer quan- 
tum Hall effect Thouless numbers defined in terms of the sensitivity of Hartree-Fock eigenval- 
ues to changes in boundary conditions, have a critical value which appears to be idependent 
of the strength of electron-electron interactions despite large changes in the density of states 
at the Fermi energy. Our work adds further motivation to efforts to clarify the connection 
between Thouless numbers and conductivities, in the presence of interactions and at a finite 
magnetic field. 
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FIGURES 

FIG. 1. Participation ratios as a function of for 7 = (crosses) and 7 = 0.4 (squares). 

The Hartree-Fock results were obtained for the case of vp = 1/2 and = 200 by averaging over 
Nd = 16 disorder reahzations. 

FIG. 2. Left panel: InP^ vs. InA^,^ for E = —0.35, — 0.6r in the absence of Coulomb inter- 
actions (7 = 0). The filling factors for these energies are 0.5, 0.066 and 0.0022 respcectively. The 
slope of the lines is —1.0 for 7^ and « 0.2 for E = Q. Right panel: Same as in the left panel 
but for 7 = 0.4. Results were calculated at the filling factors 0.5,0.081 and 0.056 for vp = 1/2 
and at 0.538,0.147 and 0.019 for up = 1/5. The Hartree-Fock eigenenergies at these filling factors 
are -0.225, -0.6 and -0.8r for up = 1/2 and -0.05, -0.40 and -0.65r for up = 1/5. The square, 
circle, and triangle symbols show results foi up = 1/2 while the pluse, cross, and diamond symbols 
show results for up = 1/5. 

FIG. 3. L^DAEiE)/N^, < 5E > /V and O.Sc/ttt^ as a function of for 7 = (left panels) and 
7 = 0.4 (right panels). L^Dae{.E) is the average number of states in the interval A£' around E 
per disorder realization. Here Np) = 1412 and up = 1/2. 
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